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ABSTRACT 

We treat resonant Raman scattering via the multiphonon exciton transi- 
tions in cylindrical quantum dots with a parabolic confinement in the lat- 
eral direction and with a finite rectangular interface-barrier confinement 
in the axial direction. The optical phonons and the electron-phonon inter- 
action are considered within the multimode dielectric model. The model 
exploits both electrostatic and mechanical boundary conditions for the 
relative ionic displacement vector, as well as the phonon spatial disper- 
sion in bulk. The confined phonon modes in a quantum dot are hybrids 
of bulk-like and interface vibrations. Raman spectra, calculated using 
the multimode dielectric model, compare well with experimental data. 
Multiphonon Raman amplitudes are calculated taking into account the 
effects of non-adiabaticity, which play a crucial role in the optical spectra 
of quantum dots. Peak intensities of Raman spectra are investigated as 
a function of the confinement frequency parameter and of the height of 
a quantum dot. Non-adiabatic transitions are shown to strongly increase 
the scattering probabilities and the relative multiphonon scattering inten- 
sities with respect to the one-phonon intensity. 



I. INTRODUCTION 



Recent experiments on Raman scattering in self-assembled quantum dots |^, ^, 
^, ^, ^ have demonstrated a rich structure of peaks provided by the optical-phonon- 
assisted quantum transitions. In contrast to bulk, specific phonon modes appear in 
quantum dots. In Ref. |^, an experimental evidence is presented, that in pyramidal 
self-assembled InAs/GaAs quantum dots, specific phonon modes exist whose frequen- 
cies differ from those both of LO and of TO phonons. They are interpreted in terms of 



interface phonons localized at the apexes of pyramids. Raman peaks assigned to in- 
terface phonons are observed in Ref . P] . Raman spectra of disk-shape self-assembled 
quantum dots, measured in Refs. 0, ^, as shown in the present paper, can be also 
interpreted in terms of phonon modes specific for quantum dots. In Ref. both 
one-phonon and two-phonon Raman bands are observed in self-assembled InAs/GaAs 
quantum dots. This result confirms the conclusion of Ref. regarding an enhanced 
efficiency of the electron-phonon interaction in quantum dots. 

The dependence of this efficiency on sizes of a quantum dot is determined by 
several factors. One of them is the following. When a finite-height potential barrier 
is at the interface of a quantum dot, an enhanced separation of the electron and 
hole charges can occur due to different masses, what enhances the efficiency of the 
exciton-phonon interaction. The second factor of this enhancement in quantum dots 
is the effect of non-adiabaticity [P]. Various theoretical investigations lead to different 
conclusions on the efficiency of the exciton-phonon interaction as a function of the 
quantum-dot size, depending on a chosen model. In Ref. [@, a donor-like exciton 
model for a spherical quantum dot has been considered, where a hole is treated as a 
localized charge in the centre of a quantum dot. It is found within this model, that 
the Huang-Rhys parameter does not depend on the quantum-dot radius. On the 



contrary, in Ref. [10|, for the same model, an increase of the excitonic polaron effect 



with decreasing quantum dot size is obtained using a variational method. The Raman 
cross section was studied for weak confinement of an exciton in a CuBr quantum dot 



in Ref. |]TT|, where it was shown that the Huang- Rhys parameter increased with the 
decrease of the dot radius. The effects of non-adiabaticity of the exciton-phonon sys- 
tem in spherical CdSe and PbS quantum dots in the strong-confinement regime were 
shown in Ref. to lead to a significant rise of probabilities of phonon-assisted tran- 
sitions and to a substantial distinction of the Raman and photoluminescence spectra 
from the Franck-Condon progression. The enhancement of the exciton-phonon inter- 
action makes possible an experimental observation of multiphonon Raman scattering 
in quantum dots, what has been confirmed by recent measurements of two-phonon 
Raman bands in self-organized InAs/GaAs quantum dots |Q. 

In the present work, we investigate multiphonon resonant Raman scattering via 
the exciton transitions in cylindrical quantum dots, which are a model for self- 
assembled disk-shape quantum dots studied experimentally in Refs. p, ^ ^. Within 
our model, we assume a parabolic confinement in the lateral direction and a finite 
rectangular interface-barrier confinement in the axial direction. This model allows us 
to calculate exciton states in a quantum dot by a variational method for a wide range 
of the lateral-confinement frequency parameter from the weak-confinement regime to 
the strong-confinement regime. 

We interpret the experiments [0, ^ ^ taking into account the effects of non- 
adiabaticity, which have been shown in Ref. to play a crucial role in the optical 
spectra of quantum dots. We demonstrate, that different additional channels of 
the phonon-assisted optical transitions open due to non-adiabaticity of the exciton- 
phonon system in quantum dots. 



II. EXCITONS AND OPTICAL PHONONS IN QUANTUM DOTS 



We consider a cylindrical quantum dot of the height h, which occupies the space 
—h/2 < z < h/2 along the 2;-axis. The lateral confinement for electrons and holes is 
ensured by a parabolic potential with the frequency parameter Qq- A medium inside 
(outside) the quantum dot is indicated by the index 1 (2). 

The exciton-phonon Hamiltonian of the system under consideration has the form 

H = Hex (re, Th) + ^Jt^ + {Puciu + Ptal) , (1) 

U V 

where H^x (fe, r/i) is the "bare" exciton Hamiltonian in the effective- mass approach, 
Te and r/i are the coordinates of an electron and of a hole, respectively; 
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with the Coulomb electron-hole interaction potential Uc (rg, rh) and with the interface 
barrier potential: 

0, \z\ < h/2, 

rh 



In Eq. (^, the notation me{h) {z) means that we take into account a distinction 
between electron (hole) effective masses in the quantum dot (medium 1) and in a 
host substance (medium 2): 



-e(.) (z) = y^'^ : ; : o' (4) 



me{h)u \z\ < h/2, 
me(h)2, \z\ > h/2. 

In the present work, we consider the case when the height of the quantum dot 
is much less than both the exciton Bohr radius and effective lateral radius of the 
quantum dot Rq ~ ^^h/ (2/iii7o), where /xi is the reduced mass of an electron and 
of a hole in the quantum dot. Eigenstates of the exciton Hamiltonian (|^), which are 
used for the calculation of the Raman spectra, are found by a variational method, 
assuming a strong-confinement regime along the z-axis. We choose the variational 
wave functions in the form 

'^NeNh,nimi, 1121712 {'^e, ^^h) = <Pn, {^e) (pN^ i^h) 

X (R) -^e'^^^^^n2m2 (P) ^e^"^^^% (5) 
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tion along the z-axis: 



where 0^ (z) are the eigenf unctions of the Hamiltonian of the electron (hole) mo- 



with the current-conserving boundary conditions at the interface 
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The functions ^mmi 9.iadipn2m2 (p) describe, respectively, quan- 

tum states of the "in-plane" center-of-mass and relative motions of an electron and of 
a hole, mi and m2 are quantum numbers of the angular momentum, while rii and n2 
are "radial" quantum numbers. The wave functions ^mmi and ipn2m2 (p) contain 
variational parameters. 

In Eq. (|l|), cUj, are frequencies of the phonon modes, which are specific for the 
quantum dot and numbered by the index z/. Pi, = 7j, (re) — •ju (ffe) is the exciton- 
phonon interaction amplitude. Within the present investigation, optical phonons in 
the quantum dot are considered using the multimode continuum model [|14|, In 
the case when the sizes of the quantum dot substantially exceed the lattice constant 
Oo, optical phonons both in the quantum dot and in the host medium are described 
by the relative ionic displacement vector u (r) = u+ (r) — u_ (r) , where u-t (r) are 
displacements of the positive and negative ions, respectively. The dynamics of the 
dispersive optical phonons in the long-wave region is determined by the generalized 



Born-Huang equation (cf. [Tg, |T7|)- It has the following form in the Fourier repre- 
sentation with respect to time: 

(^^o - u = y-^^^) E - V ■ 5 (u) , (8) 

where p is the reduced ionic mass density, E is the macroscopic electric field, uj-^ 
is the bulk TO phonon frequency at the Brillouin zone centre, ^oo and Eq are the 
high-frequency and static dielectric constants, respectively. The stress tensor 

[S (u)],^. ^ K - 24) (V ■ u) 5., + 4 (1^ + (9) 

corresponds the LO and TO bulk phonon dispersion characterized by parameters 
and fT, respectively: 



^lo (?) = ^LO (0) - vlq'' 
^TO (?) = ^To (0) - f 



(10) 



where q is the modulus of the phonon wave vector. Eigenfrequencies and basis vectors 
of phonon modes are found by the joint solution of the equation (H) and of the static 
Maxwell equation for the displacement D =£00 E + wto {eo - e^oo) pu, 

divD = 0. (11) 



For the quantum dot structure, Eqs. (||) and (|TT]) should be obeyed for each 
medium separately. The solution of this set of equations needs boundary conditions. 



They can be subdivided into electrostatic [for E (r) and D (r)] and mechanical [for 
u (r)] boundary conditions. 

Eigenvalues ixP' of the resulting boundary problem should be real. This condition 
leads to the equation, which matches the values of the functional G (u, u') = pu'-F (u) 
on both sides from the boundary between media 1 and 2, 

(G'i-G2)|,_,,,, = 0. (12) 

Here F (u) = S (u) ■ n is the force flux across the boundary, n is the unit vector 
normal with respect to the boundary. 

For structures where frequencies of optical phonons in adjacent media substan- 
tially differ from each other, optical vibrations cannot propagate from one medium to 
another. Hence, in this case we can set G\y^^^^^^^^ = 0. This approximation is equiv- 



alent to model boundary conditions proposed in Refs. ||16| , p!7| , [Tq] . The electrostatic 
interaction between lattice cells is modelled by the interaction between dipoles, which 
is strongly anisotropic. Taking into account this anisotropy, the following boundary 
conditions |T^ are applicable within this model: 

■"-L I boundary = 0' "^ll I boundary ^ ^' ^^^^ 

where subscripts _L and || denote normal and tangential (with respect to the bound- 
ary) components of a vector. For long-wavelength vibrations, the second equation of 
( |T3|) transforms into du^\ / d^nl^^^^^^^^.^ = 0, where ^„ is the coordinate directed along 
n. 

Eqs. (H) and (|ll]) are solved by substituting into them u (r) as a sum over basis 
vectors, which satisfy the boundary conditions ([T3|). As distinct from the dielectric 
continuum model (see, e. g. Ref. ^^), we find that LO phonon modes are not 



bulk-like or interface but have a hybrid character. 

From Eqs. (H) and ([TT|) with (0), dispersion equations are derived for the eigen- 
frequencies of "symmetric" (s) and "antisymmetric" (a) modes, whose electrostatic 
potentials are, respectively, symmetric and antisymmetric with respect to the reflec- 
tion in the xy-pla.iae: 

Els (?||, 1^) tanh + 62 (uj) = for symmetric modes, 

Eia {<1\\,^^) coth + 62 (to) = for antisymmetric modes, (14) 

where ( = q\\h, q\\ is the modulus of the "in-plane" phonon wave vector qy, and 62 (cu) 
is the dielectric function of a host medium. The functions eis(a) {ih^) have the form 



^ij (cih^) = ^1 (00) 
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The coefficients Xj,n are determined by 
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The summation over n in Eq. (|T3p is hmited by the value Uq ~ /i/ (2ao). This 
hmitation expresses the fact that the wavelength of an optical phonon cannot be 
smaller than twice the lattice constant oq. 

Using a formal analogy of the dispersion equations (|1^ with those for interface 
phonon modes of the dielectric continuum model (see, e. g., Ref. [0), we can 
interpret eij {q\\,uj) as the effective dielectric function of the quantum dot. As dis- 
tinct from the dielectric function in bulk, there is a set of dielectric functions in 
the quantum dot, which describe polarization response due to symmetric (j = s) 
and antisymmetric (j = a) optical-phonon modes, respectively. It is worth noting, 
that the functions eij {q\\,i^) are dispersive, since they depend both on cu and on the 
"in-plane" wave vector qy. 



III. MULTIPHONON RAMAN INTENSITIES 



Within the context of the long-wavelength approximation, the interaction of an 
electron with an electromagnetic field is described by the operator V (t) = V/e"'^^* + 
Y^^iUst^ where the terms Vj and Vg correspond, respectively, to the absorption of a 
photon with the frequency fli (incoming light) and to the emission of a photon with 
the frequency Qs (scattered light). The interaction amplitude Vk^s) is proportional 
to the projection of the electron dipole momentum operator d on the polarization 
vector e^('^) of the relevant wave: d^^^^ = e^^^^ ■ d. From the second-order perturbation 
theory, the transition probability between an initial \i) state and a final |/) state is 
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Sicofi-nj + ns); (18) 



Here, Umi and ujfi are transition frequencies. For phonon- assisted Raman scatter- 
ing in the quantum dot, \ f), and \m) are quantum states of the exciton-phonon 
system. Both the initial and the final states contain no charge carriers (electrons or 
holes), so that these states are described by a direct product of a wave function of free 
phonons with that of the exciton vacuum. Intermediate states \m) are eigenstates 
of the total Hamiltonian (^ of the exciton-phonon system. For definite polariza- 
tions of the incoming and the scattered light, the scattering probability is obtained 
by averaging Eq. (|T^ over the initial states and by summing over the final ones. 



Because electron-phonon coupling constants in such semiconductors as InAs, GaAs, 
AlAs (which are often used materials for self-assembled quantum dots) are small, the 
i^-phonon scattering intensity, corresponding to a definite combinatorial frequency 

K 

^ujy., can be analyzed to the lowest {K-th) order in the electron-phonon coupling 

i=i 

constants ai and a2- The scattering intensity is then expressed through a squared 
modulus of the scattering amphtude: 

k=l 

(19) 

Here u)^ is the frequency and rf^'''^'* = {iJi\d^^^^\Q) is the dipole matrix element of a 
transition from the exciton vacuum to the eigenstate of the Hamiltonian Hex-, f ^ 
is the inverse lifetime of an exciton in the state is the inverse lifetime of a 

phonon of the mode. 



IV. RESULTS AND DISCUSSION 



In Figs. 1 to 3, we have plotted one- and two-phonon Raman spectra for cylindrical 
GaAs/AlAs quantum dots with different values of the dot height h and of the lateral- 
confinement frequency VLq. In Fig. 1, the values of these parameters are chosen the 
same as typical values given in Ref. 0: h = 5 nm, hQ^ = 3 meV. The material 
parameters of semiconductors are taken from Ref. pT|. In particular, LO and TO 



phonon frequencies at the Brillouin zone centre are: uj^^^^ = 292.37 cm ^, u^j^^^"^ = 
268.5 cm~^, uj^^^^ = 403.7 cm~^, cj^'^'^'' = 361.7 cm~^. These frequencies are 
indicated in the figures by arrows. Raman intensities are plotted as a function of the 
Stokes frequency shift {flj — fls)- 

In order to analyze a role of non-adiabaticity, Raman spectra have been calculated 
by two methods: (i) within the non-adiabatic approach, using the Raman amplitude 
(|T9|), where transitions through all intermediate states are taken into account, and 
(ii) in the adiabatic approximation (see Ref. [^]), when non-diagonal matrix elements 
of the exciton-phonon interaction in Eq. (|19[) are neglected. The results of non- 
adiabatic and adiabatic approaches are shown in Figs. 1 to 3 by solid and dashed 
lines, respectively. 

In the one-phonon spectra, we can clearly distinguish between two groups of peaks, 
whose frequencies lie in the regions of optical-phonon frequencies of GaAs and of AlAs, 
respectively. Consequently, these peaks can be assigned to the so-called "GaAs-like" 
and "AlAs-like" phonon modes. The structure of the GaAs-like group of peaks reflects 
the spectrum of hybrid phonon modes. The hybridization of the AlAs-like phonons 
appears to be negligibly small (as distinct from the GaAs-like phonons), so that they 
are subdivided into the half-space and interface AlAs-like phonons. The interaction 
of an exciton with the half-space AlAs-like phonons is comparatively weak, because 



it is determined only by "tails" of the electron and hole wave functions outside the 
quantum dot. 

Owing to appreciable optical-phonon dispersion in GaAs, hybrid GaAs-like modes 
of the quantum dot occupy a wider frequency range than the gap between LO and 
TO phonon frequencies at the Brillouin zone centre. It appears, that a hybrid mode 
always exists, which frequency is close to the TO-phonon frequency of GaAs. This 
mode has a larger dispersion than that for other hybrid GaAs-like modes and provides 
a Raman peak of a considerable relative intensity. In the limit of large when the 
dielectric continuum model of optical phonons is applicable, the aforesaid mode turns 
into the GaAs-like interface mode. Consequently, within our model, this mode is 
formed mainly by interface vibrations and can be referred to as the "interface-like" 
mode. It is worth noting, that the experimentally observed Raman spectra of disk- 
shape self-assembled Gao.8lno.2As/Gao.5Alo.5As quantum dots [0, || (see the inset in 
Fig. 1) reveal a structure similar to that of the "GaAs-band" calculated within the 
present model. Namely, there is a comparatively narrow intense peak near uj^q^^^ 
and a wider peak between ujI^q^^^ and oo^^^\ Within our model, the former peak 
corresponds to the most-long-wavelength hybrid mode (whose amplitude is mainly 
confined to the quantum dot), while the latter is assigned to the aforesaid interface- 
like mode. 

In Fig. 2, we show Raman spectra for the quantum dot of the same height as in 
Fig. 1 but with a stronger lateral confinement: ^f^o = 10 meV. When strengthening 
the lateral confinement, the relative contributions of the interface-like phonons into 
Raman spectra diminish. In Fig. 3, we consider the quantum dot with a smaller 
height than that in Fig. 1 {h = 2 nm) and with the same value of the lateral- 
confinement frequency HQq = 3 meV. In this case, the relative intensity of peaks, 
assigned to the AlAs-hke interface modes, substantially increases. We can conclude, 
that the relative contribution of the interface-like phonons to Raman spectra is an 
increasing function of the aspect ratio R/h. 

The two-phonon Raman band consists of peaks related to different combinatorial 
optical-phonon frequencies. Note that peaks assigned to combinations of two different 
phonon modes can be higher than those related to double frequencies of each phonon 
mode (see, for example, the middle two-phonon band in Fig. 3). This fact is explained 
as follows. The light scattering by two different phonon modes can be realized through 
two channels, which differ from each other by the permutation between the phonon 
modes. Hence, the factor 2 appears for the corresponding scattering intensity when 
compared with that for the subsequent scattering by one and the same phonon mode. 

We see from Figs. 1 to 3 (where all spectra are plotted using one and the same 
scale), that the effects of non-adiabaticity considerably enhance both the absolute 
values of Raman peak intensities and the relative intensities of the two-phonon peaks 
with respect to those of the one-phonon peaks. In Fig. 4, we have plotted the 
parameter 5* = 2/2//1, (where Ii and I2 are integral intensities of the one-phonon 
and two-phonon Raman bands, respectively) as a function of the lateral-confinement 
frequency for different values of the height of the quantum dot. The parameter 5* 
is a measure of the efficiency of the exciton-phonon interaction in the quantum dot. 



Within the adiabatic approximation, S coincides with the Huang-Rhys parameter 



^^^^ = E ^ (20) 



where |/io) is the ground state of an exciton in the quantum dot. One can see, 
that S, calculated accounting for non-adiabatic transitions, strongly increases with 
strengthening confinement, and achieves the value S ~ 0.05 for hQo = 10 meV and 
h = 2 nm. The ratio of the "non-adiabatic" and "adiabatic" parameters S^^^^ / S^^\ 
as seen from the inset in Fig. 4, diminishes with increasing the confinement frequency. 
This can be explained by an increasing energy difference between various quantum 
levels of an exciton with increasing Qq, what lowers the relative probabilities of non- 
adiabatic transitions. Nevertheless, in the range of Qq studied, 5*^^"^^ > 5**^^^. 

Comparing Figs. 2 and 3 between each other, we see that with strengthening 
confinement there occurs a substantial rise of scattering intensities, especially of two- 
phonon ones. This result suggests that an experimental observation of multiphonon 
Raman bands is facilitated in small-size self-assembled quantum dots. 



V. CONCLUSIONS 



We have treated multiphonon Raman spectra in cylindrical quantum dots, tak- 
ing into account non-adiabatic effects and specific optical-phonon modes, which are 
considered within the multimode dielectric model [|T^, |r^. 

The relative contributions of different phonon modes into Raman spectra depend 
both on the quantum-dot volume and on the aspect ratio of the radius to the height 
of the quantum dot. The relative contribution of the interface- like phonons to Raman 
spectra rises when increasing R/h. 

We have shown that non-adiabatic transitions drastically enhance the efficiency 
of the exciton-phonon interaction in quantum dots. Another mechanism, which is 
responsible for this enhancement, is a separation of the electron and hole charges 
owing to different masses of an electron and of a hole, when a finite potential barrier 
exists at the quantum-dot interface. When strengthening confinement, the efficiency 
of the exciton-phonon interaction rises. 

In view of the results as obtained here, it would be important to experimentally 
investigate multiphonon Raman spectra in disk-shape self-assembled quantum dots. 
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Figure captions 



Fig. 1. One-phonon (upper panel) and two-phonon (lower panel) Raman spectra 
for a cylindrical GaAs/AlAs quantum dot with the height h = 5 nm and with the 
lateral-confinement frequency parameter hflo = 3 meV. Solid curves show results 
obtained taking into account effects of non-adiabaticity. Dashed curves show results 
of the adiabatic approximation. 

Inset: measured Raman spectra of a disk-shape self-assembled Gao.8lno.2As quantum 
dot (from Refs. Q) for different excitation frequencies. For clearness, different 
Raman spectra are shifted in the horizontal plane. The dot-dashed line indicates the 
position of the frequency shift Qj — Qs = 300 cm~^. 

Fig. 2. One-phonon (upper panel) and two-phonon (lower panel) Raman spectra 
for a cylindrical GaAs/AlAs quantum dot with the height h = 5 nm and with the 
lateral-confinement frequency parameter HQq = 10 meV. Denotations are the same 
as those in Fig. 1. 

Fig. 3. One-phonon (upper panel) and two-phonon (lower panel) Raman spectra 
for a cylindrical GaAs/AlAs quantum dot with the height h = 2 nm and with the 
lateral-confinement frequency parameter hQo = 3 meV. Denotations are the same as 
those in Fig. 1. 

Fig. 4. The parameter 5* determined as = 2J2//1, where Ii and I2 are integral 
intensities of the one-phonon and two-phonon Raman bands, respectively, as a func- 
tion of the confinement frequency Qq. In the adiabatic approximation, S is just the 
Huang-Rhys parameter. The upper group of curves shows results obtained taking 
into account non-adiabatic transitions. The lower group of curves corresponds to the 
adiabatic approximation. 

Inset: the ratio of the "non-adiabatic" and "adiabatic" parameters S^^^'^ / S^^'^ as a 
function of the confinement frequency Qq. 
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